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ABSTRACT 

The exquisite angular resolution available with Chandra should allow precision measurements of 
faint diffuse emission surrounding bright sources, such as the X-ray scattering halos created by in- 
terstellar dust. However, the ACIS CCDs suffer from pileup when observing bright sources, and this 
creates difficulties when trying to extract the scattered halo near the source. The initial study of the 
X-ray halo around GX13+1 using only the ACIS-I detector done by Smith, Edgar & Shafer (2002) 
suffered from a lack of sensitivity within 50" of the source, limiting what conclusions could be drawn. 

To address this problem, observations of GX13+1 were obtained with the Chandra HRC-I and 
simultaneously with the RXTE PCA. Combined with the existing ACIS-I data, this allowed measure- 
ments of the X-ray halo between 2-1000” . After considering a range of dust models, each assumed to 
be smoothly distributed with or without a dense cloud along the line of sight, the results show that 
there is no evidence in this data for a dense cloud near the source, as suggested by Xiang et al. (2005). 
In addition, although no model leads to formally acceptable results, the Weingartner k Draine (2001) 
and all but one of the composite grain models from Zubko, Dwek k Arendt (2004) give particularly 
poor fits. 

Subject headings: 


I. INTRODUCTION 

Practically every band of the electromagnetic spectrum 
affects or is affected by interstellar (IS) dust grains. In 
the IR, PAHs emit lines and small grains emit continuum 
radiation; in the UV/optical, small grains both extinct 
and scatter light. In X-rays, large dust grains (> 0.1 pm) 
scatter X-rays, creating halos around point sources. The 
classic paper by Mathis, Rumpl k Nordsieck (1977, 
MRN77) used the observed optical extinction to deter- 
mine the size distribution of dust grains between 0.005- 
0.25pm. Newer models, such as Weingartner & Draine 
(2001, WD01), have extended the modeling to include 
polycyclic aromatic hydrocarbons (PAHs) to match the 
observed IR emission as well as other constraints on grain 
abundances. Recently, Zubko, Dwek k Arendt (2004, 
ZDA04) found that a wide range of dust compositions 
and size distributions could fit the existing data, and 
suggested that new observational constraints from X-ray 
halos are needed to select amongst these models. 

X-ray dust scattering halos are created by the small- 
angle scattering of X-rays as they pass through dust 
grains. When an incoming X-ray interacts with the elec- 
trons in a grain large compared to the X-ray wavelength, 
the resulting Rayleigh scattering adds coherently in the 
forward direction leading to small-angle scattering; see 
van de Hulst (1957) and Mathis k Lee (1991) for de- 
tails, and Draine (2003) for a comprehensive review. 
More generally, the scattering problem can be posed as 
that of a wave interacting with a sphere, in which case 
the Mie solution applies (e.g. Smith k Dwek 1998). In 
either approach, the scattering depends largely on the 
grain size distribution, with lesser dependencies on the 
grain composition and position along the line of sight. 

Observations of X-ray scattered halos have just be- 
gun to significantly impact dust models. Smith, Edgar 
k Shafer (2002, (SES02)) described Chandra observa- 


tions of GX13+1 with the ACIS-I detector and showed 
that dust grains do not have large (> 0.8) vacuum frac- 
tions considered by Mathis k Whiffen (1989). SES02 
also found that the extremely large grains found by 
Ulysses in the solar neighborhood (Landgraf et al. 2000; 
Witt, Smith k Dwek 2001) do not seem to be common 
throughout the Galaxy. Despite these successes, SES02 
could not distinguish between the MRN77 and WD01 
models. This was in part due to calibration uncertain- 
ties as well as inherent limitations of the data. Despite 
Chandra’s excellent angular resolution, ACIS-I observa- 
tions of GX13+1 could not measure the halo within 50", 
due to massive pileup in the ACIS-I detectors. Draine 
(2003) and Xiang et al. (2005) have both noted that this 
result is therefore insensitive to dust near the source, as 
scattering from dust within the last 25% of the distance 
would lead to features primarily within the excluded 50". 
To address this shortcoming, I obtained a short Chandra 
HRC-I observation of GX13+1. The multichannel plate 
design of the HRC-I is far less sensitive to large count 
rates, which allows GX13+rs radial profile to be mea- 
sured to within 2” of the source, far closer than previ- 
ously possible. 

2. OBSERVATIONS 

GX13+1 was observed simultaneously with the Chan- 
dra HRC-I and RXTE Proportional Counter Array 
(PCA) on February 8, 2005 for 9.1 ksec (ObsID 6093) 
and 6.2 ksec (P90173), respectively. CIAO v3.3 software 
was used to process the Chandra data, which showed 
significant background flares in addition to the flux from 
the bright source. Standard processing was used for the 
RXTE data. 

2.1. Selecting Good Events 

The full-field lightcurve included significant periods 
when the count rate approached the 184 cts s -1 teleme- 



try limit. The expected HRC-I background rate for the 
full field is ~ 50cts s -11 . Despite the brightness of the 
source, the telemetry saturation was in fact primarily 
due to the particle background. After excluding a 2’ ra- 
dius circle around GX13+1, the average count rate was 

< 50 cts s -1 , but with excursions above 100 cts s -1 where 
telemetry saturation would affect the data. To eliminate 
this problem only time periods where the total counts in 
the field (i.e, > 2' from GX13+1) were < 45 cts s _1 were 
included. Although this reduced the total good time to 
3.58 ksec, ~ 200,000 counts were detected within 2' of 
GX13+1 for a source count rate of 54.8 cts s -1 . Within 
2" of GX13+1 the count rate was 42.7 cts s _1 . Accord- 
ing to §4.2.3. 1 of the Chandra Proposer’s Observatory 
Guide, the encircled energy within 2" is « 90%, rising 
to « 95% within 10". Based only on these values, it ap- 
pears that ~ 78% of the total counts are “on-axis” while 
22% are scattered by a combination of interstellar dust 
and the Chandra mirrors. 

The extremely high flux from the source combined with 
the desire to get the highest possible spatial resolution 
required an unusual instrument configuration. In col- 
laboration with the CXC Operations team, the HRC-I 
detector was positioned so that the source would appear 
in one corner of the HRC-I, while still being on-axis to 
the HRMA. This offset retained Chandra's spatial res- 
olution but ensured the source was far away from the 
normal aimpoint. Figure 1 [Left] shows the full field of 
the HRC-I, with GX13+1 at one corner. 

In Figure 1 [Right], a “jet” extending to the NE and 
containing ~ 1000 counts can be seen .•• This jet is a well- 
known detector artifact (Murray et al. 2000) which is 
normally removed by the standard processing to a level of 

< 0.1% of the total source flux 1 2 . In the case of GX13+1, 
this jet is ~ 0.5% of the apparent source count rate. The 
most likely cause is the high source count rate interfering 
with the the on-board electronic event processing (Dr. 
Michael Juda, private communication). Although the jet 
could be eliminated with aggressive filtering, this would 
also invalidate the standard calibration. It was therefore 
decided to simply ignore all events from the “jet”-side of 
the source, as shown by the box region in Figure 1 [Left] . 

2.2. Extracting the Spectrum and Flux 

The surface brightness of the X-ray halo must be nor- 
malized by the source flux to make absolute measure- 
ments of the dust column density. GX13+1 is observed 
almost constantly since it is a RXTE All-Sky Monitor 
(ASM) source with a average rate of 20-30 cts/s. How- 
ever, as the RXTE ASM has little spectral sensitivity 
and the HRC-I has no effective energy resolution, simul- 
taneous RXTE PCA observations of GX13+1 were taken 
to obtain a useful spectrum of the source. Although the 
RXTE PCA itself has only moderate resolution and lit- 
tle sensitivity below 2 keV, GX13+Ts spectrum is dom- 
inated by 2-4 keV photons (SES02). The PCA spectrum 
is shown in Figure 2 [Left] as fit with a simple model 
consisting of an absorbed multi-color disk model plus 
a blackbody, following Ueda et al. (2004). The col- 
umn density was fixed at the value found by Ueda et al. 
(2004) from the Chandra HETG, H = 3.2 x 10 22 cm- 2 , 

1 http://asc.harvard.edu/proposer/POG/html/HRC.html 

2 http://cxc.harvard.edu/cal/Hrc/Documents/ghosts.ps 


since this result is far more accurate than one obtained 
from the PCA. The best-fit inner temperature of the mul- 
ticolor disk and the blackbody were 1.73 keV and 3.52 
keV, with absorbed 1-10 keV fluxes of 7.43 x 10~ 9 erg 
cm _2 s _1 and 2.30 x 10 _1 °erg cm _2 s -1 , respectively. De- 
spite the large reduced y 2 (> 50, driven by systematic 
errors), this model is an adequate fit for this work since 
only the total flux and the approximate spectral shape 
are needed to calculate the predicted response of the 
ChandraHRC-I. Nonetheless, it is important to note that 
fluxes measured with the RXTE PCA are systematically 
high by 10-15% in the 2-10 keV range, compared to other 
X-ray observatories 3 . 

To check the expected count rate, we folded this spec- 
trum through the HRC effective area file for on-axis 
Cycle 7 data (hrciD2005-ll-30pimmsN0008.fits), as 
shown in Figure 2 [Right]. The total predicted source 
count rate in the HRC-I is 64.6 cts s -1 , ~ 18% larger 
than the observed HRC-I count rate of 54.8 cts s _1 within 
2’ of GX13+1. The discrepancy is primarily due to the 
overestimation from the RXTE PCA calibration, with an 
additional complication due to spatial variation in the re- 
sponse of the HRC-I that reduces the effective area of the 
detector corners relative to the center 4 . 

The Chandra PSF, measured as a ratio of the sur- 
face brightness to the source flux, is the background for 
this observation. The RXTE PCA, a non-imaging detec- 
tor, includes both the direct source flux and the scat- 
tered halo photons, which must be removed to avoid 
double-counting. However, as the goal is to measure 
the scattered halo fraction itself, this problem is recur- 
sive. I addressed this by assuming a column density of 
3.2 x 10 22 cm -2 and calculating the total scattered frac- 
tion for the MRN77, WD01, and ZDA04 BARE-GR-B 
models, weighted by the HRC-I response. The resulting 
halo fraction ranged from 13-26%. This predicted halo 
strength is consistent with the result that 22% of the to- 
tal source counts are between 2" — 120". Therefore, for 
purposes of calculating the background PSF, the RXTE 
PCA flux was reduced by 20% to exclude the halo contri- 
bution, with a 7% systematic error. This reduction is in 
addition to the 15% reduction described above. The 7% 
error is likely not the dominant term in the systematic er- 
ror, however. The observation of a very bright source in 
one corner of the HRC-I detector is at the extreme edge 
of the available calibration, and so careful consideration 
of all uncertainties will be required. 

Another concern regarding this observation was that 
a significant short-term change in the source flux, on 
the order of 12-24 hours, would also affect the halo in 
a time-delayed manner ( e.g . Vaughan et al. 2004). At 
smaller angles the delay could be even longer. The RXTE 
ASM data was checked for a 10 day period before the 
observation, but no strong or significant variation was 
seen. Although Type I X-ray bursts have been seen from 
GX13+1 which show 3 — 4x the normal flux, they only 
last ~ 15 seconds (Matsuba et al. 1995). In this case, no 
bursts were seen in either the HRC-I or PCA lightcurves, 
and indeed the halo observation would not be sensitive 
to such a small variation. 

3 http://universe.gsfc.nasa.gov/xrays/programs/rxte/pca/flux.scale. 

4 See http://asc.harvard.edu/cal/Hrc/flatfield.html 
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Fig. 1. [Left] The full HRC-I observation of GX13+1, with the excluded rectangular region marked. [Right] Expanded image of 
GX13+1 on the HRC-I, showing detector “jet” and excluded region. The 2” radius circle shows the near-source region excluded due to 
likely detector non-linearity. 
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Fig. 2. [Left] The RXTE PCA spectrum of GX13+1, fit .with an 
model folded through the HRC-I response. The integrated count rate 

2.3. Point-spread function 

An accurate measurement of the Chandra HRC-I 
point-spread-function (PSF) between 2-100” from the 
source is crucial to this observation. An accurate ray- 
trace model (ChaRT 5 ), of the Chandra HRMA has 
been calibrated for near-source (< 2") photons , SES02 
showed that at large scattering angles this model sig- 
nificantly underpredicts the PSF, leading to substantial 
problems in the analysis. Therefore, SES02 relied upon 
an ACIS-I observation of Her X-l as a PSF calibrator, 
but this source is affected by pileup within ~ 10” and 
therefore cannot be used in the 2-10” range. 

As shown in Figure 2 [Right], the spectrum of GX13+1 
peaks at ~ 2 keV. The HRC-Fs lack of spectral response 
means that spectral differences between any calibration 
source and GX13+1 will lead to additional complica- 
tions. The best possible calibration source would be 
a bright, hard, and lightly-absorbed X-ray source ob- 
served on-axis with the HRC-I. The X-ray binary LMC 
X-l matches these requirements reasonably well, and two 
Chandra observations (ObsID 1200, 1201) of the source 

J http://asc.harvard.edu/soft/ChaRT/cgi-bin/www-saosac.cgi 



absorbed multicolor disk plus blackbody model. [Right] The best-fit 
predicted is 64.6 cts/s. 

have been done. However, they were both done early in 
the mission (August 1999) before the HRMA final focus 
was set and are thus unsuitable. Since then, the bright- 
est hard X-ray source with little absorption and a known 
(albeit variable) flux to be observed with the HRC-I is 
3C273 (ObsID 461 on Jan 22, 2000). Figure 3 shows 
3C273’s surface brightness, divided by its source flux, as 
observed with the Chandra HRC-I (excluding the well- 
known jet region) . The spectrum was taken from a Chan- 
dra HETG observation done twelve days earlier (ObsID 
459) which is well-fit by an absorbed power-law with T = 
1.67±0.01 and Fy(2 — lOkeV) = (1.08±0.03)xlCT 10 ergs 
cm~ 2 s -1 . The absorption column was fixed at the Galac- 
tic value, Nh = 1-8 x 10 2O cm -2 . The predicted HRC- 
I count rate (based on the CXC PIMMS tool) for this 
spectrum is 8.7 cts/s, while the actual source count rate 
was 26% higher at 11 cts/s. As the source is vari- 
able, this was taken as showing little change and the 
flux was simply assumed to have increased by 26% dur- 
ing the HRC-I observation. The core of the PSF was 
fit with a Gaussian term centered at 0 with FWHM of 
1. 007^0 005” and amplitude 1676 ± 42 arcinin -2 . In ad- 
dition, the best-fit model included two power-law terms 
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Fig. 3. [Top] The radial profile of 3C273’s surface brightness, 
divided by the source flux, fit with the sum of a Gaussian plus two 
power laws and a constant. [Bottom] The ratio of the data and 
model, showing small excursions at 1.5" and 10", but generally 
good agreement. 


TABLE 1 

Smooth Dust Model Parameters 


Model 

Nh(HRC) 
10 22 cm~ 2 

Nh(ACIS) 
10 22 cm -2 

xl 

MRN77 

2.4 ±0.2 

2.85 ±0.05 

2.0 

WD01 

1.51 ±0.02 

2.0 ±0.2 

3.6 

BARE-GR-S 

2.7 ±0.2 

2.9 ±0.1 

1.9 

BARE-GR-FG 

2.6 ±0.2 

3.00 ± 0.04 

1.9 

BARE-GR-B 

3.6 ±0.3 

3.4 ±0.4 

2.8 

BARE-AC-S 

2.5 ±0.2 

3.0 ± 0.1 

2.1 

BARE-AC-FG 

2.5 ±0.1 

3.0 ±0.2 

2.2 

BARE-AC-B 

3.3 ± 0.3 

3.60 ± 0.05 

2.2 

COMP-GR-S 

2.02 ± 0.02 

2.9 ±0.4 

4.9 

COMP-GR-FG 

2.23 ±0.04 

3.0 ±0.3 

3.5 

COMP-GR-B 

3.0 ±0.2 

3.71 ±0.09 

1.9 

COMP-AC-S 

2.41 ± 0.02 

3.7 ±0.6 

7.0 

COMP-AC-FG 

2.67 ± 0.02 

3.8 ±0.5 

5.0 

COMP-AC-B 

4.24 ± 0.04 

6.7 ±0.9 

5.7 

COMP-NC-S 

11.1 ± 0.6 

13.5 ±0.9 

2.4 

COMP-NC-FG 

2.81 ±0.02 

4.6 ±0.8 

8.6 

COMP-NC-B 

3.39 ± 0.03 

5.9 ± 1.1 

11.7 


with Ti = 4.06 ± 0.05, r 2 = 2.401 q o 2 an d amplitudes 

Ai — (3.4toJ) x 10~ 7 , A 2 = (2.12to.o!) x 10 -4 arcmin -2 
at 100”. The particle and sky background was fit with 
a constant, (2.60 ± 0.02) x 10~ 3 arcmin -2 . As Figure 3 
shows the fit is quite good over a large range of surface 
brightnesses, with the somewhat large reduced xt = 2.6 
likely due to the extreme precision of the measurement 
compared to the relatively simple model. 

3. RESULTS 

The HRC-I observations are most useful between 2- . 
100” , since beyond that radius the ACIS-I data can 
measure the energy-resolved X-ray halo. Therefore, the 
ACIS-I data were reprocessed (with CIAO 3.3) and rean- 
alyzed following the approach described in SES02 except 
as noted below. Both the HRC-I and ACIS-I results were 
used in the final analysis. We note that in reprocessing 
the ACIS data, the source flux measurement, done via 
the CCD transfer “streak” , was redone with a better cal- 
ibration and improved handling of the background sub- 
traction which resulted in an overall ~ 15% decrease in 
the measured source flux. The calibration changes in- 
clude a spatially-varying modification of order ±5% in 
the quantum efficiency uniformity in CALDB 2.28, and 
an energy-dependent increase of up to 16% in the overall 
effective area which was added in CALDB v3.2.1. 

The data were fit using the CIAO fitting engine Sherpa 
using scattering models based on the exact Rayleigh- 
Gans (RG) approximation (Smith & Dwek 1998). 
This model assumes the grains are spherical but uses 
the energy-dependent optical constants rather than the 
Drude approximation when calculating the scattering ef- 
ficiency. Smith & Dwek (1998) noted that the full Mie 
treatment is necessary for X-rays < 2keV, since the 
RG overestimates the total scattering at low energies. 
GX13+l's spectrum, however, is dominated by photons 
with E > 2keV, even for the predicted spectrum ob- 
served by the HRC-I (see Figure 2 [Right]). Therefore, 
the simpler RG treatment is justified. 

The initial analysis assumed the dust was smoothly 
distributed along the line of sight. Unlike SES02, where 
the predicted PSF was subtracted from the data, here the 
PSF was incorporated into the fitting directly to allow 


for an explicit inclusion of uncertainty in the PSF. Fits 
to the ACIS-I and HRC-I data included a constant fac- 
tor that allowed for calibration uncertainty in the overall 
PSF, caused primarily by systematic errors in the source 
flux. For both the ACIS-I and HRC-I data this multiplier 
was allowed to vary by up to 10%. In many cases the fit 
pushed the multiplier to an extremum of the range, show- 
ing that systematic uncertainties remain in the data, al- 
though it is not clear what component dominates them. 
In SES02, systematic errors in the PSF manifested as 
energy-dependent column density fits, since to first or- 
der an error in the PSF could be adjusted by changing 
the overall halo scattering. As the total halo intensity 
is inversely proportional to energy, this effect is often 
a linear dependence of the best-fit Nh on energy. To 
check for this, I allowed the value of Nh to vary inde- 
pendently in the HRC-I and each energy band of the 
ACIS-I data. I used the F-test to determine that in only 
one case (ZDA04 BARE-GR-B) were the best-fit ACIS-I 
Nh values better described by a linear energy-dependent 
model than by a constant value. Even in this case, the 
F-test significance was only 3.5%, a negligible value given 
the number of different models tried. It seems unlikely, 
therefore, that there is a significant error in the relative 
power in the dust-scattered (halo) and mirror-scattered 
(PSF) photons. 

Table 1 shows the best-fit Nh results for the HRC-I and 
the “average” ACIS-I value fit and the total xl assuming 
smoothly-distributed dust along the line of sight for the 
MRN77, WD01, and the 15 ZDA04 models. These can be 
compared to the value of 3.2 x 10 22 cm -2 found by Ueda 
et al. (2004). The ACIS-I column densities are ~ 20% 
larger than the HRC-I values, although the models with 
the lowest \l values tend to have the best agreement. 
The most likely cause of this discrepancy is cumulative 
errors in the source flux measurements combined with 
calibration differences between the ACIS-I and HRC-I 
detectors. 

Figure 4 shows the profile of the HRC-I data along 
with the ACIS-I data at 2.5 ± 0.1 keV, near the median 
energy of the spectrum as observed by the HRC-I. This 
figure shows the level of agreement between the HRC-I 
and ACIS-I data agree with each other in the overlap re- 
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gion (50 — 100"), as well as the large (> 20x) difference in 
the HRC-I and ACIS-I backgrounds in the 500 — 1000" re- 
gion. The radial profile shown in Figure 4 is fit assuming 
the line of sight (LOS) dust is “smoothly-distributed” 
and has a composition and size distribution described 
by WD01 [Left], and the ZDA04 BARE-GR-S [Right] 
models. Both models generally agree with the data, al- 
though the WD01 model undepredicts the ACIS-I data 
in the 300 — 500" range, while conversely the BARE-GR- 
S model underpredicts the HRC-I data in the 10 — 50" 
range. 

The smoothly distributed dust fits show small dis- 
crepancies that might be indications of dusty molecular 
clouds along the line of sight. These would appear as 
“bumps” in the profile whose position and strength de- 
pends upon the relative distance to the cloud and its col- 
umn density. I therefore refit the HRC-I and ACIS-I data 
using a two-component model that included a smoothly- 
distributed component plus a cloud with variable posi- 
tion and column density. It is worth noting that mod- 
els using only a single cloud with no smooth component 
gave generally poor fits independent of the dust model 
used, and so were not considered in detail. This is not 
unexpected since (a) SES02 was unable to fit a single 
cloud using only the ACIS-I data and (b) GX13+1 is 
reasonably near the Galactic center ((Lb) — (13.5°, 0.1°), 
D = 7 ± lkpc; Bandyopadhyay et al. (1999)) where a 
sightline dominated by a single cloud would be unusual. 

To reduce fit time, the column density of both the 
smooth component and a cloud was fixed to be the same 
for all datasets, as was the position of the cloud along the 
line of sight. While more realistic than allowing the cloud 
column density to vary as a function of X-ray energy, this 
has the effect of magnifying residual systematic errors. 
As noted previously, the halo strength diminishes with 
energy while the relative PSF strength increases which 
can create a trend in the best-fit column density as a 
function of energy. However, since in only one case out 
of fourteen was such a trend seen previously, it seems un- 
likely that the systematic errors are driving the resulting 
best-fit parameters in the smooth plus cloud model. 

The best-fit parameters for each model are shown in 
Table 2. None of the fits are formally acceptable (x 2 
ranges from 2.1 to 9.9)', although some are clearly better 
than others. In cases where the best-fit position is 0, the 
la upper limit is shown. 

4. DISCUSSION 

Smith, Edgar & Shafer (2002) analyzed the ACIS-I 
observations of GX13+Ts dust scattered halo and found 
that the dust size distribution does not extend to very 
large (> l^m) grains and that grains do not have a large 
vacuum fraction. However, the ACIS-I data could not 
distinguish between the MRN77 and WD01 models, as 
the two distributions lead to similar scattering profiles at 
large angles. Similarly, the data left open the possibility 
that there might be a substantial population of grains 
near the source (Draine 2003; Xiang et al. 2005). These 
would create a near-source scattered halo that was ob- 
scured by pileup in the ACIS-I detector. The primary 
goal of the HRC-I observation was to remove these un- 
certainties by measuring the halo near the source. This 
would determine which dust model best fit the data, as 
well as detecting (or put limits upon) variation in the 


TABLE 2 

Smooth Plus Cloud Model Parameters 


Model 

Nh (smooth) 
10 22 cm -2 

Nn(cloud) 
10 22 cm -2 

Relative 

Position 

xi 

MRN77 

2 60^ u ' ua 
z - ou -0.05 

0 69 +u ‘ ut * 
u.oy_ 0 03 

< 0.003 

2.4 

WD01 

1.38 ± 0.02 

0.46 ± 0.01 

< 0.001 

3.5 

BARE-GR-S 

3.02 ± 0.01 

0.15 ±0.01 

0.89 ± 0.01 

2.1 

BARE-GR-'FG 

2.95 ±0.01 

0.11 ±0.01 

0.91 ±0.01 

2.2 

BARE-GR-B 

3.51 ±0.01 

0.57 ±0.01 

0.82 

2.4 

BARE-AC-S 

2.6 ±0.1 

0.21 ±0.06 

< 0.001 

2.4 

BARE-AC-FG 

2.6 ±0.1 

0.23 ±0.06 

< 0.001 

2.4 

BARE-AC-B 

3.4 ±0.1 

0.35 ±0.14 

0.85 

2.1 

COMP-GR-S 

1.6 ±0.1 

0.82tg;{2 

< 0.0003 

2.9 

COMP-GR-FG 

2.0 ±0.1 

0-62t°;“ 

< 0.004 

2.6 

COMP-GR-B 

3.2 ±0.1 

0.23 ± 0.05 

< 0.01 

2.2 

COMP-AC-S 

1 £+ 02 

1.3 ±0.1 

< 0.0002 

3.7 

COMP-AC-FG 

o O + 0- 2 
*'*-0.1 

1.0± 

< 0. 

3.3 

COMP-AC-B 

2.9 ±0.2 

2.3 ±0.1 

< 0.0004 

3.6 

COMP-NC-S 

11 4 +0 - 6 
AA '^-0.3 

1.3 ±0.3 

jO.0008 

2.5 

COMP-NC-FG 

1 6 +0>3 

1.9 ±0.2 

< 0.0001 

4.3 

COMP-NC-B 

1 6 +0<9 

2 e i +0 - 4 

< 0.00001 

6.5 


dust distribution along the line of sight. 

The fit results shown in Tables 1 and 2 contain a few 
surprises. Just as in Smith, Edgar & Shafer (2002), 
the WD01 model had the smallest column density of any 
of the models when fit with either smoothly distributed 
dust or after adding a cloud. However, the overall re- 
sult was a significantly worse fit than found with either 
the MRN77 or many of the ZDA04 models. As Figure 4 
shows, the smooth WD01 model fits the HRC data well, 
but underpredicts the halo measured by ACIS between 
150" — 400", while the MRN77 model underestimates the 
halo measured by the HRC between 10 — 50". Examining 
the other models show that these two cases are represen- 
tative. Adding a single dust cloud to the model results 
in a solution with a cloud 70-90% of the distance to the 
source if the pure smooth model underestimates the halo 
around 30". Conversely, adding a cloud component to 
smooth dust models that underestimate the halo around 
300" tend to put the cloud near the Sun. 

Although uncertainties remain due to calibration un- 
certainties, the lack of energy resolution in the HRC-I, 
and inability of X-ray data alone to distinguish between 
dust models, the overall quality of the fits shown in Ta- 
ble 1 and Figure 4 do not support the proposition that 
a significant cloud of dust is present near GX13+1. This 
result is confirmed by the fits shown in Table 2, which 
suggest that if a cloud is present it is either very near 
the Sun or 70 - 90% of the distance to GX13+1 with 
Nh ~ 5 x 10 21 cm -2 . In any event, the cloud contributes 
less than half of the total LOS column density as mea- 
sured solely via the X-ray halo. This result therefore dis- 
agrees with results of Xiang et al. (2005), who analyzed 
the radial profile of the HETG observation of GX13+1 
and found that more than 70% of the dust was effectively 
“at" the source using either the MRN77 or WD01 grain 
models. This type of distribution would appear as a large 
increase in the surface brightness between 2 — 10", which 
is simply not apparent in our data. 

ZDA04 described in detail how constraining dust mod- 
els requires combining multiwavelength data from the 
IR to X-rays while simultaneously considering the metal 
abundances in the grains. Due to the nature of op- 
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Fig. 4. [Left] The X-ray halo from smoothly-distributed WDOl-type grains fit to the radial profile of GX13-|-l’s surface brightness 
divided by the source flux. The HRC-I observations are fit with a red line; the ACIS-I data (at 2.5 keVj with a blue line. [Right] Same, for 
the ZDA04 BARE-GR-S model. 


tical/UV extinction and X-ray scattering, few sources 
show strong signatures of dust in all of these wave- 
bands (Valencic & Smith 2007). Nonetheless, it is 
possible to constrain the allowed dust models by com- 
paring the column density predicted by the models to 
that measured using other techniques. In the case of 
GX13+1, measurements of the column density range 
from 2.5 — 4.0 x 10 22 cm -2 (Charles & Naylor 1992). 
Optical measurements provide only an upper limit of 
2.9 x 10 22 cm -2 based on plausible but unconfirmed 
assumptions about the source spectrum(Garcia et al. 
1992) . The total H I column density through the Galaxy 
at the position of GX13+1 is 1.8 x 10 22 cm -2 (Dickey & 
Lockman 1990) , but this is misses the contribution from 
molecular H 2 that is likely to be strong in the Galac- 
tic plane. The HETG observation of GX13+1 agrees 
(weakly) with these results (Ueda et al. 2004), al- 
though it does not strongly limit it. Only Mg can be 
directly measured ( Nm 9 = l-Bdl^g x 10 18 cm -2 ), equiv- 
alent to Nh= 4.8ti'^ x 10 22 cm -2 assuming solar abun- 
dances. The 2 a upper limits for Si and S are equivalent 
to Nh < 4 x 10 22 cm -2 . Of course, LMXBs have shown 
significant variable internal absorption (Hertz & Grind- 
lay 1983) in X-rays, so this spectral measurement sets at 
best an upper limit to the actual interstellar component 
that is responsible for the halo. 

Despite these difficulties in independently measuring 
the total LOS dust column density, we can reasonably 
justify excluding the value of Nh = (1.11 ± 0.06) x 
10 23 cm -2 found in Table 1 for the ZDA04 COMP-NC- 
S model fit to the HRC-I data. However, this was only 
model of theirs that used composite grains without bare 
carbon grains that had a plausible value of ■ Although 
more data are needed, this class of models, along with 
the group of “composite grains with bare amorphous car- 
bon" models are clearly suspect since they do not gener- 
ate an X-ray halo similar to these observations. In fact, 
the only smoothly-distributed ZDA04 composite grains 
model that fit with y 2 < 3 had graphitic carbon and 
B star abundances (COMP-GR-B). After adding a dust 
cloud to the model, the COMP-GR-B model fit with 
X 2 = 2.2 while the next best fit (excluding the unrealistic 
COMP-NC-S model) was COMP-GR-FG with = 2.6, 
a significantly worse fit. 


5. CONCLUSIONS 

The principal results from this analysis are: 

1. Although challenging, HRC-I observations can be 
used to recover the near-source region excluded by 
pileup in the ACIS-I detector. The lack of energy 
resolution can be finessed if another measurement 
of the source spectrum is available. 

2. Fitting the source profile and background PSF in- 
dependently improves overall results, since calibra- 
tion uncertainties in the flux from the source and 
background objects can then be included explicitly. 

3. There is no sign of a significant near-source dust 
cloud in the radial profile, as suggested by Xiang 
et al. (2005). 

4. In agreement with SES02, the WD01 model under- 
predicts the total column density to the source, and 
again leads to poor fits, although not so bad as to 
be excluded given the calibration uncertainties. 

5. Some models from the ZDA04 paper, if not con- 
clusively excluded, are at least implausible. In 
general, the ZDA04 models with composite grains 
(excepting the graphitic carbon model with B star 
abundances) gave poor fits, while the bare carbon 
and silicate grain models tended to fit well. 

It should be noted that the relatively good fits found 
using the simple smoothly-distributed dust model are 
somewhat surprising, since X-ray halos probe both the 
largest grains whose size and composition are the least 
constrained from observations in other wavelengths. Ad- 
ditionally, X-ray halos are primarily observed through 
highly-absorbed lines of sight. These probe dust in dense 
molecular clouds that cannot be observed in the optical 
or UV due to the extremely large extinction. Finally, 
all of these models assume spherical grains, although re- 
cently some calculations have been done on nonspherical 
grains (Draine & Allaf-Akbari 2006) that show the ef- 
fects are small (< 20%) for the WD01 model at ~ 2keV. 
Despite these potential problems, a number of existing 
grain models agree quite well with the observations, sug- 
gesting grains in dense clouds (with the exception of the 
densest regions that take up very little volume and may 
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be optically-thick to X-rays) are not too dissimilar from 
grains in less dense regions. 

I thank Michael Juda for his substantial assitance in 
arranging the observation and understanding the HRC- 


I calibration. I would also like to thank Terry Gaetz 
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Chandra PSF. Finally, I also thank Eli Dwek for many 
helpful discussions. This work was supported by the 
NASA Chandra observation grant G05-6144X. 
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